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Short time approximation is developed for system-environmental bath mode interactions involving 
a general non-Hermitian system operator A, and its conjugate, A*, in order to evaluate onset of 
decoherence at low temperatures in quantum systems interacting with environment. The developed 
approach is complementary to Markovian approximations and appropriate for evaluation of quantum 
computing schemes. Example of a spin system coupled to a bosonic heat bath via A oc <r_ is worked 
out in detail. 
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I. INTRODUCTION 

Quantum system exposed to environmental modes is 
described by the reduced density matrix, and its evolu- 
tion deviates from the ideal, usually pure-state, dynam- 
ics. For short times, appropriate for quantum computing 
gate functions and, generally, for controlled quantum dy- 
namics, approximation schemes for the density matrix 
have been developed recently |lH3j- The present work 
derives a new rather general short-time approximation 
which applies for models with system-bath interactions 
involving a general system operator. It thus extends the 
previously known approach 00 which was limited to 
couplings involving a single Hermitian system operator. 

We consider an open quantum system with the Hamil- 
tonian 



H = H S + H E 



(1) 



Here Hs describes the system proper. It is coupled to 
the environment (bath), described by Hb- The system 
and bath are coupled by the interaction Hi. The bath 
has been traditionally modelled by a large number 
of uncoupled bosonic modes, namely harmonic oscillators 
(with ground state energy shifted to zero), 



H B = 



(2) 



Here bk are the bosonic annihilation operators corre- 
sponding to the bath modes, and from now on we use 
the convention H = 1. 

In most of this work, we consider the general system 
Hs, and we assume that the interaction with the bath 
involves the system operator A that couples linearly [T3 - 
IbH to the bath modes, 



H I =Aj2g k bl+A^9kh, 



(3) 
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with the interaction constants gk- 

Let R(t) denote the overall density matrix. It is com- 
monly assumed |llll4j that at time t = the system and 
bath are not entangled, and the bath modes are thermal- 
ized, 



R(o) = P (o)l[e k 



where 



ik = 



z i e _^ fe6 t fc 



with j3 = 1/kT, and 



(4) 



(5) 



(6) 



We point out that while the quantum system S, de- 
scribed by the reduced density matrix pit), is small, typ- 
ically two-state (qubit) or several-qubit, the bath has 
many degrees of freedom. The combined effect of the 
bath modes on the system can be large even if each of 
them is influenced little by the system. This has been 
the basis for the arguments for the harmonic approxi- 
mation for the bath modes and the linearity of the 
interaction, as well as for the Markovian approximations 
[Tfjlil^ l that assume that the bath modes are "reset" to 
the thermal state by the "rest of the universe" on time 
scales shorter than any dynamical time of the system in- 
teracting with the bath. 

The frequencies of the oscillators of the bath are usu- 
ally assumed to be distributed from zero to some cutoff 
value u> c . The bath modes with the frequencies close 
to the energy gaps of the system, A-Ey — Ef — Ej , con- 
tribute to the "resonant" thermalization and decoherence 
processes. Within the Markovian schemes, the diagonal 
elements of the reduced density matrix of the system, 



pit) =Tr flj R(t) 



(7) 



approach the thermal values oc e~ Ei l kT for large times 
exponentially, on time scale T\. The off-diagonal el- 
ements vanish, which represents decoherence, on time 
scale T2, which, for resonant processes, is given by T2 — 
1T\. However, generally decoherence is expected to be 
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faster than thermalization because, in addition to reso- 
nant processes, it can involve virtual processes that do 
not conserve energy. It has been argued that this ad- 
ditional "pure" decoherence is dominated by the bath 
modes with near-zero frequencies |T(1T3T^ . At low tem- 
peratures, this "pure decoherence" is expected 0] to 
make T 2 < T x , 

Since the resetting of these low-frequency modes to 
the thermal state occurs on time scales 1/kT = (3, the 
Markovian approach cannot be used at low tempera- 
tures p!ollllll5| . Specifically, for quantum computing in 
solid-state semiconductor-heterostructure architectures 
[I6H23I ]. temperatures as low as few 10 mK are needed. 
This brings the thermal time scale to f3 ~ 10~ 9 sec, which 



is close to the single-qubit control times 10 



-li 
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16-23]. Alternatives to the Markovian approximation 
have been suggested pill^ . 

In this work, we generalize the recently suggested 
scheme , applicable for Hermitian A only, to a wider 
class of interaction Hamiltonians. We treat the case when 
the system operator A entering the interaction, see is 
not Hermitian. In actual applications in quantum com- 
puting, calculations with only a single qubit or few qubits 
are necessary for evaluation of the local "noise," to use 
the criteria for quantum error correction |30H35| . For ex- 
ample, the system Hamiltonian is frequently taken pro- 
portional to the Pauli matrix a z . The interaction opera- 
tor A can be proportional to a x , which is Hermitian. Such 
cases are covered by the short-time approximation devel- 
oped earlier 00- However, one can also consider models 
with A cx <7_. Similarly, models with non-Hermitian A 
are encountered in Quantum Optics |ll| . In Section II, 
we develop our short time approximation scheme. Re- 
sults for a spin-boson type model are given in Section 
III. 



II. SHORT TIME APPROXIMATION 

In this section we obtain a general expression for the 
time evolution operator of the system lll.il) within the 
short time approximation. The system operators H$ and 
A need not be specified at this stage; the derivation is 
quite general. 

In order to define "short time," we consider dimension- 
less combinations involving the time variable t. There 
are several time scales in the problem. These include 
the inverse of the cutoff frequency of the bath modes, 
l/w c , the thermal time (3 — 1/kT, and the internal char- 
acteristic times of the system 1/AEij. Also, there are 
time scales associated with the system-bath interaction- 
generated thermalization and decoherence, Ti^- The 
shortest time scale at low temperatures (when (3 is large) 
is typically 1/uj c . The most straightforward expansion in 
t yields a series in powers of u> c t. The aim of d eve loping 
more sophisticated short-time approximations |2l3| is to 
preserve unitarity and obtain expressions approximately 
valid up to intermediate times, of order of the system 



and interaction-generated time scales. The latter prop- 
erty can only be argued for heuristically in most cases, 
and checked by model calculations. 

The overall density matrix, assuming time- 
independent Hamiltonian over the quantum-computation 
gate function time intervals fl6ll23l ] , evolves according to 



R(t) = U(t)R(0) [U(t)?, 



where 



U(t) 



-i{H s + H B +H,)t 



(8) 



(9) 



is the evolution operator. 

The general idea of our approach is the following. We 
break the exponential operator in © into products of 
simpler exponentials. This involves an approximation, 
but allows us to replace system operators by their eigen- 
values, when spectral representations are used, and then 
calculate the trace of R(t) over the bath modes, obtain- 
ing explicit expressions for the elements of the reduced 
density matrix of the system. For Hermitian coupling op- 
erators, A^ = A, our approach reduces to known results 



We split the exponential evolution operator into terms 
that do not have any noncommuting system operators in 
them. This requires an approximation. For short times, 
we start by using the factorization 36-38] 



e ~i(H s +H B +H,)t+0{t 3 ) 
e -%H s t e -i(Hi+H B )t e -iH s t 



(10) 



where we have neglected terms of the third and higher 
orders in t, in the exponent. The middle exponential in 

035, 



-i{H B +Hi)t _ e -i(H B +AG t +A t G)t 



where 



(12) 



still involves noncommuting terms as long as A is non- 
Hermitian. In terms of the Hermitian operators 



L S i(A + At), 
M=i(A-At), 



(13) 
(14) 



we have 

AG t + A t G = L(G + G t ) +iAf (G-G*). (15) 

We then carry out two additional short-time factoriza- 
tions within the same quadratic-in-i (in the exponent) 
order of approximation, 

3 = e i[M(G-&)-iH B ]t e t,H B t ( 16 ) 
x e -i[H B +L(G+G- i )]t e ±H B t e ^[M(G-G r )-iH B }t_ 
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This factorization is chosen in such a way that H remains 
unitary, and for M = or L = the expression is iden- 
tical to that used for the Hermitian case 00 • The evo- 
lution operator then takes the form 



U 



- c-i.Hst- 



■iH s t 



with 5 from 1(16(1 . which is an approximation in terms of 
a product of several unitary operators. 

It has been recognized 2.3] that approximations of this 
sort are superior to the straightforward expansion in pow- 
ers of t (or more exactly, u c t). Specifically, in (|1U|) . we 
notice that Hg is factored out in such a way that Hb, 
which commutes with H$, droppes out of all the commu- 
tators that enter the higher-order correction terms. This 
suggests that a redefinition of the energies of the modes 
of Hb should have only a limited effect on the corrections 
and serves as a heuristic argument for the approximation 
being valid beyond the shortest time scale l/w c , up to 
intermediate time scales. 

Our goal is to approximate the reduced density ma- 
trix of the system. We consider its energy-basis matrix 
elements, 



where 



Pmn(t)=TY B (m\UR(0) rf\n) 



H s In) = E n In) 



(18) 



(19) 



We next use the factorization 1)10116(1 to systematically 
replace system operators by c-numbers, by inserting de- 
compositions of the unit operator in the bases defined by 
Hs, L, and M. First, we collect the expressions 1(41161171 
Hl^l , and use two energy-basis decompositions of unity to 
get 

pq 



X 



Tr B [(m\E\p)l[e k (q\Et\n)]. (20) 



operators to c-numbers. For convenience, let us define 
the operators 



7T Q/37 = | a) (a\ (3) (0 \i) (7I 



(23) 



(17) and 



U S2 S3 = e s l9'k b k t + 8 29 k '>lt-is 3 u, k blb k t^ / 2 4) 



The resulting expression for the trace entering l(20|) is 

1t B {(m\E\p)l[e k (q\E^n)} (25) 

k 

fj,j Xj A.' 

where the indices A and p label the eigenstates of L and 
M, respectively, and 

r k = Tr k ^Ui tll _i tll iU 00 _iU- lXl - l \ 1 ,i (26) 



1 o,o, 



"2 2^2,-2^2,2 



1^3, 2^3,- 3 



X ^0,0 i^iA2,iA 2 ,-lWo,0,i^-i M 4,|M4,-|} • 

In order to calculate the trace over the fcth bath mode 
in ((2*f]|l . we rearrange the operators using the cyclic prop- 
erty, in such a way that the formula ((Al(l . derived in 
Appendix A, can be used to simplify products of two or 
three operators Wata time. For example, we can transfer 
the operator Ui fllt _t fll> 1 to the right hand side, getting 
the combination 



1,,. _iWi,„ _i„, i- 

2 , 2 "4 1 2 2 "^ ' 2 r* 1 ' 2 



(27) 



inside the trace, and use the identity I(A1(I . We then 
transfer U 00 _i to the right hand side and repeat the 
process, now for the three rightmost operators. After 
several steps we arrive to the following expression for the 
trace, 



We now define the eigenstates of L and M, 

L|A) = A|A), (21) where 

M |m) = H |M) ■ (22) 

u = 

The operators 3 and E< introduce exponentials in l(20|) 
that contain either L or M in the power. By appro- 
priately inserting \ty (M or 1^) (a*I between these 

A M 

exponentials, we can convert all the remaining system and 



7fc — {0kU u ,-u*,o} ■ 



(28) 



i ( _^ki \ „ 2i _*^k± . uj k t 
e 2 — l #t" e 2 Sln \_ 

LO k t\ ) p ~ 



i 

UJ k t 



UJ k t 



e 2 l e 2 _ 1 j 



(29) 
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TZk = exp • 



i Is* 



,sm- —J + + (smLj k t - uj k t) \-\ + 



4 sin — sin — /U_/x + -2sm — (A_/x + - A+^_) j j 



(30) 



Here we introduced the variables 

M2 ± M3: 



and 



A± — Ai ± A2 



(31) 
(32) 



(33) 



The trace in l|28l) can be evaluated, for instance, by using 
the coherent states technique, see Appendix iBl 



*e 



<-\g k \ 2 t 2 coth^ 



The expression which follows from Q25[l. I|29|) . (|30|l and 
(Ell is 



(34) 

where 



H% = exp[-V (*)], 



(35) 



V = B 2 (t) (\_ + nLfi'L) + B 2 (t/2) (m'I - M'-) - ^ (*) (M- - M-) A_ 

- »C (t) A_A+ - iC (t/2) + (jfifi'l) + iS (t) (A_a»+ - A +i u'_) - iCi (t) //_//| 

I 



(36) 



The coefficients here are the spectral sums over the bath 
modes, 



B 2 {t) = 2 y\9k\ sin 2^i coth ^ 

; ,,,2 O O 



2~ 



C(t) = 2%i-(u-A/-Miia,v./i: 



* w * 



(37) 



(38) 



these functions are well known |39||40j. The result also 
involves the new spectral functions 



S(*) = -2£^Lsin^, 



(39) 



F (*) = 4 Z^ _ T 2- sm ^- sm ^- coth ^-- (40) 



* w * 



4 " 2 2 

Furthermore, for the sake of convenience we defined 

Ci (t) = 2C (t/2) - C (t) . (41) 

By using l|2U|) and (|35|l . we obtain our final result for 
the density matrix evolution, 

Pm»(t) = ^^ei^+V^-y VM(0 ) (42) 
x (m| 7r^ lAl/i2 |p) (g| 7r AI3 A 2M4 |n) e~ p , 



where the first sum over p and g is over the energy eigen- 
states of the system; the second sum is over Ai, A2 and 
Hi, . . . , fJ-4, which label the eigenstates of the operators L 
and M, respectively; see (|2"H) and (|2"2|) . 



III. DISCUSSION AND APPLICATION 



The result (|42|l looks formidable in the general case. 
However, in most applications evaluation of decoherence 
will require short-time expressions for the reduced den- 
sity matrix of a single qubit. Few- and multi-qubit sys- 
tems will have to be treated by utilizing additive quan- 
tities |4lll43j . accounting for quantum error correction 
(requiring measurement), etc. For a two-state system — a 
qubit — the summation in l|42|) involves 2 8 = 64 terms, 
each a product of several factors calculation of which is 
straightforward. Still, the required bookkeeping is cum- 
bersome, and we utilized the symbolic language Math- 
cmatica to carry out the calculation for an illustrative 
example. 

We consider the model 0] defined by 



(43) 
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where A > is a constant, a± — \{cr x ± &y) and tr z 
are the Pauli matrices, b k and b k are the bosonic cre- 
ation and annihilation operators, and gu are the cou- 
pling constants. Physically this model may describe, for 
example, a qubit interacting with a bath of phonons, or 
a two-level molecule in an electromagnetic field. In the 
latter case, this is a var iant of the multi-mode Jaynes- 
Cummings model [ll|45j- Certain spectral properties of 
this model, the field-theoretic counterpart of which is 
known as the Lee field theory, are known analytically, 
e.g., 0| . However, the trace over the bosonic modes, to 
obtain the reduced density matrix for the spin, has not 
been obtained exactly. 



For the model 143(1 we have A = er_ and = a + , 
so that L — <r x /2 and M = a y /2. We have | A12) = 
(|t) ± \[)) /v2j with eigenvalues Ai,2 = ±1/2, and 
iMi.a) = (IT) ±* ID) A/2, with eigenvalues Hi t % = ±1/2. 
For the initial state, let us assume that the spin at t = 
is in the excited state |t) (T|, so that the initial density 
matrix has the form 

P(0)=(JS). (44) 

Calculation in Mathematica yields the following results 
for the density matrix elements, p\2{t) = and 



4pu (0 = 2+ e - 2B2{t) + e- iB2 (^ cosh (2 F) + 2 e - 2B2 ^) sinh (B\) cos (5) (45) 
+ 2e- B2 (5) C os(Ci)sin(5) + ie- B2(t) - B2 (^) [e iCl sinh (—iS + F) + e~ iCl sinh (—iS — F)] , 

I 



where C\ was defined in 141|) and 

B x (t) = 2B 2 {t/2) - B 2 (t) . (46) 

Where not explicitly shown, the argument of all the spec- 
tral functions entering ((43|l is t. 

In order to obtain irreversible behavior and evaluate a 
measure of decoherence, we consider the continuum limit 
of infinitely many bath modes. We introduce the density 
of the bosonic bath states V (ui), incorporating a large- 
frequency cutoff u) c , and replac e the sum mations in P7)l- 
lUS} by integrations over u> [^ Ildl39l47| . For instance, 
(p?7j) takes the form, 

B* { t)= L.^M^sin^coth^. (47) 
J uo z 2 2 

o 

We will use the standard Ohmic-dissipation |j expres- 
sion, with an exponential cutoff, for an illustrative calcu- 
lation, 

V(u) \g(cu)\ 2 = (48) 

where Q is a constant. 

We point out that the results obtained for the den- 
sity matrix elements depend on the dimcnsionless vari- 
able w c t, as well as on the dimensionless parameters f2 
and L0 C P (= HuJc/kT, where we remind the reader that H, 
set to 1, must be restored in the final results). Interest- 
ingly, the results do not depend explicitly on the energy 
gap parameter A, see (|43|) . This illustrates the point 
that short-time approximations do not capture the "res- 
onant" relaxation processes, but rather only account for 
"virtual" relaxation/decoherence processes dominated by 
the low-frequency bath modes. However, the short-time 
approximations of the type considered here are mean- 
ingful only for systems with well-defined separation of 




FIG. 1: Schematic behavior of s (t) for different values of Q, 
decreasing from i to iv. 

the resonant vs. virtual decoherence processes, i.e., for 
H/A ^ l/u c . For such systems, fij A = l/A defines one 
of the "intermediate" time scales beyond which the ap- 
proximation cannot be trusted. 

As an example, we calculated a measure of deviation of 
a qubit from a pure state in terms of the "linear entropy" 

S (f)=l-Tr [ p 2(t)]. (49) 
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FIG. 2: The comparison between the 0(t 2 ) expansion, i, and 
the short-time approximation, ii. 



Figure 1 schematically illustrates the behavior of s(t) for 
different fl values, for the case ui^ 1 << (3. The val- 
ues of s(t) increase from zero, corresponding to a pure 
state, to 1/2, corresponding to a completely mixed state, 
with superimposed oscillations. For Ohmic dissipation, 
three time regimes can be identified |4jj. The short- 
est time scale is set by t < 0(l/u} c ). The quantum- 
fluctuation dominated regime corresponds to O (l/w c ) < 
t < 0(l/kT). The thermal-fluctuation dominated 
regime is t > O (1/kT). Our short time approximation 
yields reasonable results in the first two regimes. For 
t > 0(l/kT) it cannot correctly reproduce the process 
of thermalization. Instead, it predicts approach to the 
maximally mixed state. 

Figure 2 corresponds to the parameter values typical 
for low temperatures and appropriate for quantum com- 
puting applications, lo c (3 — 10 3 , with f2 = 1.5 • 10~ 7 
chosen to represent weak enough coupling to the bath 
to have the decoherence measure reach the threshold for 
fault-tolerance, of order 10 -6 , for "gate" times well ex- 
ceeding 1/oJc, here for ui c t over 10. The leading-order 
quadratic expansion in powers of the time variable t is 
also shown. Its validity is limited to t < O(l/oj c ) and 
it cannot be used for evaluation of quantum-computing 
models. 

This research was supported by the National Security 



Agency and Advanced Research and Development Activ- 
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APPENDIX A 

Our aim is to derive a relation of the form 

where the operator U si ,s 2 ,s 3 was defined in Consider 
the quantity 

A = e K ( bt + a i)( b + Q 2) e 7ifc t +72fc e -2 ; (fc t +/3i)(fc+/32) ; ( A2 ) 

where , b are the bosonic creation and annihilation op- 
erators, x, oti, (3i are c-numbers. Let us use the identity 

e *b-ptff ^ 6 t) e -*b+,3bt = / (& + ft &t + a ) , (A3) 

to represent the first and third exponentials in A in the 
form 

e x(b^ +a 1 ~)(b+a 2 ) _ ^aib—aib^ e xb^ b e -aib+a 2 b i (A4) 
e -z(&+ + /3i)(&+/32) _ gPib-fctf e -xrfb e -Pib+Pib\r^\ 



We then combine the second exponential in (|A2|) and the 
last and first exponentials in i|A4|l and (|A5(1 . by utilizing 
the identity 



(A6) 



which follows from (|A3|) . The resulting exponential op- 
erator, with exponent linear in b and 6 T , is sandwiched 
between e xb b and e~ xb b . Therefore, the following iden- 
tity can be utilized [ill], 



xb+b 



^-xb 1 b 



= f{be- x ,tfe x ) 



(A7) 



Once again using (|A6f) , we arrive at the following expres- 
sion, 



where 



vb+fj,b' +r 



fi = {a 2 -fa) (e x - 1) + jie x , 
v = (fa - ai ) (e- x - 1) + l2 e~ x , 



(A8) 



(A9) 
(A10) 



and 
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r = — 2 {a.\p2 — 012P1) sinh — + (0110:2 — P1P2) sinhx 
+ ("l + A) (e x - 1) + ^72 (02 + #2) (e _a! - 1) . 



Now l|Al|l follows, with 



k = exp 



x exp 



x exp 



2 \g k \ . 2 / xu^A 

sin I — — I - W2^i ) 



x 2 w 2 



X 2 Lul 



(sin (xui k t) (viWi - v 2 w 2 ) + xuj k t (v 2 w 2 - Ui«?i)) 



2a;ijjfe 



r 



(All) 



(A12) 



and 



'7 : 



XLO k t 

i 1 



(e M - 1) («! + v 2 ) + W3 e"" fct , (A13) 



XLOkt 



-ixukt 



1) (wi+w 2 ) + u; 3 e^ fct .(Al4) 



where the integration in complex plane is defined via 

(B5) 



1 



d z — —dxdy. 

7T 



APPENDIX B 

Let us calculate the trace in (|28|l which has the form 

T = Tr je^e 1 *^ 61 '}, (Bl) 

where we omitted the index k since all the calculations 
here are in the space of a single mode. We use the 
coherent-state technique 0|. The coherent states \z) 
by definition are eigenstates of the annihilation operator, 
b 



b\z)=z\z) 



(B2) 



with complex eigenvalues z = x + iy. These states are 
not orthogonal 

(zi| z 2 ) = exp [zlm - i |zi| 2 - i \z 2 \ 2 ^ , (B3) 

and they form an overcomplete set. The identity operator 
can be written as 



d 2 z\z) (z\ = 1, 



(B4) 



We represent the trace (JBlf) by the coherent-state in- 
tegral using the relation 



TrA = J d 2 z(z\A\ 



(B6) 



where A is an arbitrary operator. We then use the normal 
ordering, Af, formula for bosonic operators, represented 
schematically (see f° r details) by 



(B7) 



The second term in the trace in (|B1() is split by using 
(|A6|) . All instances of b and can then be replaced by z 
and z* , and the integral evaluated to yield the expression 
for the trace, 



Hp- coth 4 



T 



(B8) 
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